Tunnel ionization of open-shell atoms 
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A generalized ADK (Ammosov-Delone-Krainov) theory for ionization of open shell atoms is com- 
pared to ionization experiments performed on the transition metal atoms V, Ni, Pd, Ta, and Nb. 
Our theory is found to be in good agreement for V, Ni, Pd, and Ta, whereas conventional ADK 
theory overestimates ionization by orders of magnitude. The key to understanding the observed 
ionization reduction is the angular momentum barrier. Our analysis shows that the determination 
of the angular momentum barrier in open shell atoms is nontrivial. The Stark shift is identified 
as the second dominant effect responsible for ionization suppression. Finally, these two effects can- 
not explain the Nb data. An analysis of the electron spins suggests that Pauli blocking might be 
responsible for the suppression of tunneling in Nb. 
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Tunnel ionization theory was developed in several semi- 
nal works 0, Q • All of these theories are based on the 
single active electron (SAE) approximation, where only 
the weakest bound electron interacts with the laser held. 
The most commonly used tunneling theories are ADK 
(Ammosov-Delone-Krainov) and MO-ADK (molecular- 
ADK) theory 0, which are in excellent agreement with 
experiments in noble gases and small molecules 

The situation is different for more complex systems. 
Theoretical predictions from ADK theory often overesti- 
mate experimentally measured ionization yields by orders 
of magnitude |J, |5( . Recently, a quasi-analytical multi- 
electron theory of tunnel ionization was developed Q , re- 
vealing that the SAE approximation loses its validity in 
complex materials. A comparison with ionization exper- 
iments in Cgo showed reasonable agreement. This theory 
opens the door to explore the new multi-electron physics 
associated with tunnel ionization in complex systems. 

In this paper we develop the theory of Ref. [f| fur- 
ther to analyze tunnel ionization of the transition metal 
atoms Vanadium (V), Nickel (Ni), Palladium (Pd), Tan- 
talum (Ta), and Niobium (Nb). They are significantly 
more complex than noble gas atoms due to an open d 
shell and large angular and spin momenta. ADK theory 
overestimates ionization in transition metal atoms by 
orders of magnitude. Our work reveals the following new 
findings: 

(i) We find good agreement with V, Ni, Pd, and Ta ex- 
periments. This further corroborates the applicability of 
the theoretical framework, developed here and in Ref. 
0, to general complex systems. 

(ii) The first major effect responsible for the suppression 
of ionization is found to be the angular momentum bar- 
rier. Due to the complexity of open shell atoms, determi- 
nation of I is not as straightforward as in noble gases, but 
requires a careful analysis. This reveals unprecedented 
insight into the structure and complexity of transition 
metal atoms. The second dominant contribution to the 
suppression of tunneling comes from an increase of the 



ionization potential caused by the Stark shift. Our anal- 
ysis demonstrates the importance of the Stark shift for 
tunnel ionization. 

(iii) The two effects discussed in (ii) cannot explain the 
strong ionization reduction found in Nb. An analysis of 
the spins of s- and d-orbital electrons suggests that the 
reduction of tunnel ionization in Nb might be caused by 
Pauli blocking, an effect central to many-body physics. 

(iv) Starting from the quasi-analytical multi-electron the- 
ory of of Ref. 0, a generalized multi-electron ADK/MO- 
ADK theory is derived, which includes correction factors 
accounting for multi-electron effects and angular momen- 
tum barrier. 

Our analysis begins with the tunnel ionization rate in 
atomic units, as given in Ref. 0, 
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where I and m refer to the angular momentum and mag- 
netic quantum numbers of the field- free asymptotic wave- 
function of the weakest bound valence electron that is 
matched with the wavefunction of the tunneling electron. 
In the case that m is a non-degenerate, good quantum 
number, e.g. homonuclear diatomic molecules, w = w m , 
with w the ionization rate. When m is not a good 
quantum number, several m contribute to the tunnel- 
ing wavefunction and w = J2 m Wm - ^or m degenerate, 
e.g. noble gas atoms, the ionization rate is determined by 
w = 1/(2Z + 1) J2 m Wm - Note that Eq. (JIJ corrects a mis- 
print in Eq. (8) of Ref. with respect to the position of 
Further, B hn = Y^m'^ln'mi^Qlm'Qm', and Qi m 
is the normalization constant of the spherical harmonics 
0, 0. The rotation matrix V l m , m rotates the coordi- 
nate system by the Euler angle R, and allows calcula- 
tion of ionization along an arbitrary direction, once the 
matching coefficients Ci m have been determined along 
one direction. Finally, t = J z 1 p z dz, I/A 2 = J 1 l/p z dz, 
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The electric field E is assumed to point along the z- 
direction, Z is the charge of the residual ion after ion- 
ization, and — p n and (3 + — (3 n -\ refer to the polariz- 
ability of the n-electron system before and after ioniza- 
tion [3. The Stark shifted Q and the field free ioniza- 
tion potential of the weakest bound electron are given 
by k 2 /2 =I p -(P- /3 + )E 2 /2, and I p , respectively. The 
four terms in Eq. represent the barrier coming from 
the Coulomb potential, the laser polarization, the image 
charge, and the angular momentum, respectively. The 
last three terms are not accounted for in ADK /MO- ADK 
theory. Integration is performed between the limits Zq 
and Zi, where z% is the outer turning point. 

The lower integration limit Zq is the point at which the 
laser dressed, tunneling wavefunction is matched to the 
field free, asymptotic ground state, see Ref. [||. Match- 
ing has to be done at a field free point. As complex sys- 
tems are highly polarizable, the polarization field cancels 
the laser field at a point a, where fl+E/a 2 — aE w 0. 

1/3 

Solution of this equation gives zq = a = pJ ' . Due to 
shielding the inner part of the system, z < a, is approxi- 
mately field free. In a metallic sphere, the matching point 
lies on the surface. 

When V(z) < k 2 /2 for z < z < z\ is fulfilled, the 
square root in p z0 can be Taylor expanded. In this limit, 
all the integrals can be solved analytically and contact 
can be made between Eq. JTJ) and ADK/MO-ADK the- 
ory. The calculation is fairly involved; details will be 
given in a longer paper. The generalized ADK/MO-ADK 
ionization rate is given by 
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where w refers to the conventional ADK/MO-ADK ion- 
ization rate, 
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The ADK/MO-ADK matching coefficient in Eq. © is 
related to the matching coefficient in Eq. Q by Ci m = 
Ci m e KZ » Z Q Z/K . We use here w = 1/(21 + 1) £ m w m . The 
three correction terms in a come from the angular mo- 
mentum barrier, the laser polarization barrier, and the 
image charge term, respectively. The 1-barrier is not con- 
tained in the ADK/MO-ADK theory. The ^-dependence 
in ADK/MO-ADK theory comes from the spherical har- 
monics and the matching coefficient. In the limit of I = 
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FIG. 1: Tunnel ionization yield J wdt in V, Ni, Nb, Pd, and 
Ta versus laser peak intensity. The open circles denote the 
experimental results from Ref. 0- The dashed and full lines 
refer to the theoretical results obtained from ADK theory, and 
from our theory, Eq. (JjJJ, respectively. Gaussian laser pulse 
parameters: wavelength 1.5/im and FWHM width 50/s. 



and small polarizability, Eq. (J2J) reduces to the conven- 
tional ADK/MO-ADK ionization rate, w m = w . A com- 
parison of Eq. © with Eq. shows excellent agree- 
ment for the transition metal atoms investigated here. 
Therefore, we use Eqs. © - @ throughout the paper. 
The applicability of © for complex molecules will be 
investigated in a future work. 

Key to understanding tunnel ionization in transition 
metal atoms is knowledge of the angular momentum / 
of the asymptotic wavefunction of the weakest bound va- 
lence electron. Quantum chemistry codes have been used 
to determine I for small molecules and noble gas atoms 
0, Q • As Hartree-Fock (HF) is sufficient to model these 
systems, the asymptotic wavefunction can be obtained 
from the highest occupied molecular orbital. We find 
that this method is also applicable to Pd and gives 1 = 2, 
see the discussion below. This allows us also to extract 
Ci, n for Pd from the numerical analysis and to compare 
it with the ADK matching coefficient 0, as determined 
by quantum defect theory. The ADK matching coeffi- 
cient is commonly used to calculate tunnel ionization of 
atoms. The numerical result is found to be by a factor of 
1.2 smaller. This difference is insignificant and corrobo- 
rates the validity of quantum defect theory. Hence, the 
ADK matching coefficient is used for all the other tran- 
sition metal atoms, where a numerical determination is 
not possible. 

In all other transition metals the wavefunction under- 
goes a complex configuration mixing, when the valence 
electron is pulled away from the nucleus to the asymp- 
totic region. Quantum chemistry codes are not sensitive 
to the asymptotic part of the wavefunction, as they op- 
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timize the binding energy, which mainly depends on the 
wavefunction close to the nucleus. Although an accurate 
description of the asymptotic wavefunction is in principle 
possible within a configuration interaction (CI) analysis, 
in practice it requires an unreasonably large number of 
configurations. Therefore, due to the lack of proper nu- 
merical methods I is inferred from the configuration data 
given in table [I] 

Our analysis of the transition metal atoms starts with 
Vanadium. Whereas in noble gases / is determined by the 
highest occupied orbital, this is no longer possible here. 
Although the highest occupied orbital of V is a 4s orbital, 
the following more careful analysis of the total angular 
momentum L shows that I = is forbidden. The total 
angular momentum of the atom and ion ground state is 
L = 3 and L + = 2, respectively, see table H] The plus 
subscript refers to ions and ion-parameters throughout 
the paper. From the angular momenta L and L + a range 
of possible values of I can be derived. This is done by us- 
ing the field free, asymptotic, atomic ground state, which 
can be approximated as the tensor product of the ionic 
ground state and of the valence electron far away from the 
nucleus. As a result, the total angular momentum of the 
ion plus valence electron, as determined by the angular 
momentum addition theorem, has to be equal to the to- 
tal angular momentum of the neutral atom. This is only 
possible when I fulfills the condition L + — 1<L< L + + / , 
yielding 1 < / < 5. The possible range of / can be fur- 
ther narrowed down by utilizing the fact that parity also 
has to be conserved. As V and V + have even parity, the 
parity of the asymptotic valence electron wavefunction 
must be even, too. This requires I to be even and thus, 
/ = 2,4. As the 3d orbitals are energetically close to 
the highest occupied 4s orbitals, the dominant tunneling 
channel comes from a 3d-orbital with 1 = 2. 

We would like to note three points associated with the 
above analysis. First, the conservation laws were applied 
to the field free, asymptotic ground state, which is the 
initial state for tunneling. Angular momentum and par- 
ity are of course not conserved during tunnel ionization. 
Second, the atomic ground state close to the nucleus and 
in the asymptotic limit is described by different config- 
urations. This mixing of configurations in the time in- 
dependent ground state should not be confused with a 
time-dependent reconfiguration. Third, we assume that 
the atom relaxes adiabatically during ionization into its 
ionic ground state. When an excited state is populated 
after ionization, which is termed shake-up || , the ioniza- 
tion potential increases by the energy difference between 
ion excited and ground state, AI. In transition metal 
atoms AI w 0.3 — leV is small, so that the shake-up 
channel might be important. Whereas it would be sur- 
prising when tunnel ionization leaves ions dominantly in 
an excited state, due to the complexity of the problem 
it also cannot be rigorously ruled out. Our angular mo- 
mentum analysis gives I = 0, 2 for the shake-up channel. 



In case I = 2, the shake up channel gives a small correc- 
tion to the fundamental channel and our results remain 
unchanged. In case I = 0, our theory would predict a too 
large ionization yield, close to ADK theory. The I = 
shake-up channel is excluded based on the fact that I = 2 
explains the experimentally observed ionization suppres- 
sion very well. We suggest an experiment to corroborate 
our conclusion. The first excited state in transition metal 
atoms is non-dipole allowed and therewith metastable. 
Therefore, the population of the ion states after ioniza- 
tion can be probed by exciting dipole allowed transitions 
and measuring fluorescence or the loss experienced by 
the probing laser. For example, in V+ the lowest dipole 
allowed transitions from the 3d 4 ground state and the 
first excited (metastable) 3d 3 4s state can be probed, with 
transition energies of ~ 4.5eV and 4eV, respectively. In 
dependence on whether the shakc-up channel has / = or 
2, one would find a predominant population of the 3ci 3 4s 
state or of the 3d 4 ground state, respectively. 

Given the simple approach used here and experimen- 
tal uncertainties, our theory is in good agreement with 
the experimental data of Ref. @. In Fig. [2 the ion- 
ization yield J w[E(t)]dt is plotted versus the laser peak 
intensity. For pulse parameters and shape of the laser 
electric field E(t) see the figure caption. Our theory (full 
line) explains the difference between experiment (empty 
circles) and ADK theory (dashed line). Half of the differ- 
ence is due to the /-barrier which is not accounted for in 
ADK/MO-ADK theory. As /?+ is small, the role of the 
laser polarization induced barrier is insignificant. The 
other half of the difference is accounted for by the Stark 
shift, which reduces ionization by increasing the ioniza- 
tion potential. The large Stark effect arises from the big 
difference in polarizabilities f3 and /3+, see table [I] As 
the Stark effect increases with laser intensity, its effect 
is most pronounced at high intensities. Therefore, the 
ionization yield rises gradually and reproduces the shape 
of the experimental curve, in contrast to the abrupt rise 
predicted by ADK/MO-ADK theory. Switching off the 
1-barrier shows that the Stark effect requires the 1-barrier 
to suppress ionization effectively. Without the 1-barrier 
ionization would be shifted to lower intensities, where the 
role of the Stark effect is insignificant. 

Part of the large discrepancy between ADK/MO-ADK 
theory and experiment might arise from the fact that 
ADK/MO-ADK overestimates tunnel ionization in the 
above barrier limit. The Coulomb barrier is suppressed 
at an electric field Ef, s = Ip/4, which corresponds to a 
laser intensity of 7 x 10 12 W/cm 2 in V. Hence, barrier 
suppression in V occurs before ADK/MO-ADK theory 
shows that significant ionization takes place. In H and 
He, ADK/MO-ADK theory overestimates above the bar- 
rier ionization rates by a factor of 2-4 |1Q], which is by far 
not enough to explain the difference in Fig. ^ In tran- 
sition metals the magnitude of the overestimation is un- 
known. Nevertheless, we believe that barrier suppression 
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is not the dominant factor responsible for the difference 
between ADK/MO-ADK theory and experiment. Signif- 
icant ionization usually takes place around the intensity 
at which the barrier is suppressed. In V ionization be- 
comes dominant in the range between 4 x 10 13 W/ cm 2 and 
7 x 10 13 W/cm 2 , which is by a factor of 5-10 larger than 
the barrier suppression intensity. It appears to be coun- 
terintuitive that ionization should only start at intensities 
so much higher than the barrier suppression intensity. In 
our opinion the large difference points to the fact that 
ADK/MO-ADK theory, which is exclusively based on the 
Coulomb barrier, underestimates the true tunneling bar- 
rier. This is corroborated by our theory, Eq. (J2J). Due 
to ^-barrier and Stark shift, the tunneling barrier is not 
suppressed over the laser intensity range used in Fig. 

For Ni and Pd we get similarly good agreements as 
for V, see Fig. For Ni our analysis yields I = 2. 
In Pd determination of I is straightforward. The total 
angular momentum changes from L = to L + = 2. The 
configuration changes from 4d 10 to 4d 9 , showing that the 
tunneling electron comes from a d-orbital and that 1 = 2, 
in agreement with the HF analysis discussed above. 

In Ta it was not possible to unambiguously identify 
I. From table [3 we find two relevant ionization channels. 
The possible angular momentum values are I = 0, 2, 4, 6. 
Hence, ionization from an s-orbital and from a d-orbital 
are possible. In V tunneling from a d-orbital is preferred, 
probably because breaking of the closed subshell 4s 2 con- 
figuration is energetically unfavorable. We suspect that 
tunneling from a d-orbital in Ta might be favored for 
similar reasons. This is supported by comparing theory 
and experiment for 1 = 2, which yields good agreement. 
For I = the experimental data cannot be explained. 

Finally, for Nb we could not get agreement with ex- 
perimental data. The possible values of I are I = 0, 2, 4. 
Niobium does not have a closed s-subshell like V, Ni, and 
Ta, which might favor tunneling of a d-electron. There- 
fore, we use I = in Fig. We have also performed test 
calculations for I = 2 and 4. These also do not give sat- 
isfactory agreement, corroborating our conclusion that 
the Z-barrier is not responsible for the observed ioniza- 
tion suppression. Further, ionization due to dimers and 
higher clusters was excluded experimentally [l2^ . 

A spin analysis of the s- and d-electrons shows that the 
spins of the tunneling and bound electrons are all parallel 
in Nb. In all other transition metal atoms the spins of 
the bound electrons are parallel, however the tunneling 
spin points in the opposite direction. Therefore, in V, Ni, 
Pa, and Ta, the exchange effect between tunneling and 
bound electrons is zero. In contrast to that, the tunneling 
electron in Nb feels the exchange effect with the bound 
electrons. The exchange effect makes it harder for the 
tunneling electron to pass the other open shell electrons 
and to escape. In the language of many-body physics 
this is termed Pauli blocking. 



Based on the fact that none of the other mechanisms 
can explain the Nb experiments, we speculate that Pauli 
blocking is responsible for the observed ionization sup- 
pression. A quantitative proof will require a generaliza- 
tion of the current ionization theory to account for the 
exchange effect, which will be subject to future work. 
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TABLE I: Table for the transition metal elements V, Ni, 
Pd, Ta and Nb, containing the total charge of the nucleus, 
configuration, total spin and angular momentum, polariz- 
ability (A 3 ) of neutral and singly charged ion, ionization 
potential I p and the angular momentum I of the tunneling 
electron. The plus subscript refers to singly charged ions and 
their properties throughout the paper. The polarizability 
is calculated by using density functional theory in the 
local density approximation |llf : f3 agrees with published 
data, see http:// webbook.nist.gov; for /3+ no data is avail- 
able for comparison. The configuration data is taken from 
http://physics.nist.gov/PhysRefData/Handbook/index.html 
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